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EXACT COMPUTATION OF THE n-LOOP INVARIANTS OF KNOTS 


STAVROS GAROUFALIDIS, ERIC SABO, AND SHANE SCOTT 


Abstract. The loop invariants of Dimofte-Garoufalidis is a formal power series with arith¬ 
metically interesting coefficients that conjecturally appears in the asymptotics of the Kashaev 
invariant of a knot to all orders in 1/A^. We develop methods implemented in SnapPy that 
compute the first 6 coefficients of the formal power series of a knot. We give examples that 
illustrate our method and its results. 
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1. Introduction 


1.1. The Volume Conjecture to all orders in 1/N. The best known quantum invariant 
of a knot in 3-space is the Jones polynomial [Jon87]. The Kashaev invariant {K)is[ of a 
knot K (for N = 1,2,...) [Kas95] coincides with the evaluation of the Jones polynomial 
of a knot and its parallels at complex roots of unity [MMOl]. The Volume Conjecture of 
Kashaev [Kas97] states that for a hyperbolic knot K, 


lim —log I (A) AT I = 
N^oo 7V 


Vol(A) 
27r ^ 


where Vol(A) is the hyperbolic volume of K. An extension of the Volume Conjecture to all 
orders in 1/A was proposed independently by Gukov and the hrst author [Guk05, Gar08]. 
Namely, for every hyperbolic knot K there exists a formal power series 0iy(h) G C[h] such 
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that 

( 1 ) {K)n ~ iV3/2e^^^0^(27r*/iV), 

where Ck is the complexihed volume of K divided by 27ii, 

( 2 a) (pxifi) = 

( 2 b) 0 +(h) ei + hFKlhl 

( 2 c) tk G Fk, 

and Fx is the trace held of K. 

1.2. Ideal triangulations, shapes, and the loop invariants. The left hand side of 
Equation (1) is concretely dehned given a planar projection or an ideal triangulation of a 
knot, and is typically given by a hnite state-sum where the summand is a ratio of quantum 
factorials. Examples of state-sum formulas for the Kashaev invariant of the 4i, 62 and 61 
knots are given in [Kas97, (2.2)-(2.4)]. 

On the other hand, the power series 4>k{K) that conjecturally appears in the right hand 
side is not an explicit function of the knot. Numerical computations of the Kashaev invariant 
were performed by Zagier and G., and using numerical interpolation and a variety of guessing 
methods, it was possible to recognize the hrst few coefficients of the power series (j)K for 
several knots [GZ]. 

The main result of Dimofte-G. [DG13] was the construction of a power series 0.y(h) that 
depends on an ideal triangulation of a knot complement. For a detailed discussion on ideal 
triangulations and their gluing equations, see [Thu77, NZ85] and also [DG13, Sec.1.2]. Ex¬ 
plicitly, an ideal triangulation T with N tetrahedra gives rise to a vector z = (zi,..., z^) of 
shapes that satisfy the Neumann-Zagier equations 

N 

j=i 

ioT i = 1,..., where z" = X — Xjzy These equations are obtained as follows. Fix an edge of 
the ideal triangulation T, and set the product of the shape parameters of all tetrahedra that 
go around the edge e equal to 1. If ^ is the shape of a tetrahedron that contains the hxed 
edge, then its contribution to the above product is 2 ; or z' or z" according to the convention 
of Figure 1. Finally, replace z' = —( 22 ")“^ using the relation zz'z" = —1. This gives rise to 
the equations (1.2), one for each edge of T. Likewise, there is an equation os the same type 
for each peripheral curve. 

In the above equations, we have removed one edge equation and replaced it with the 
meridian cusp equation. Neumann-Zagier [NZ85] prove that (A|B) is the upper part of a 
symplectic matrix. It follows that (A|B) has rank A and AB^ is symmetric, where B^ is 
the transpose of B. We will assume that B has nonvanishing determinant. Furthermore, 
we will assume that our triangulation T is such that there exists a solution to the gluing 
equations in (C \ {0,1})'^ that recovers the complete hyperbolic structure of the hyperbolic 
knot K. In that case, z is a vector of algebraic numbers and Q(^i,..., z^) is a number held 
(the shape held) which coincides with the invariant trace held and with the trace held of the 
knot [NR92, Thm.2.2,2.4]. Finally, one can choose a hatenning that is an integer 
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Figure 1. An ideal tetrahedron and its shape assignment. 

solution of the linear equation A/ + B/" = v. This determines a Neumann-Zagier datum 
7 = f, which in turn dehnes the power series (j)'y{h). Of course, different 

ideal triangulations give rise to different Neumann-Zagier data, hence to potentially different 
formal power series On the other hand, the left hand-side of Equation (1) depends 

only on the hyperbolic knot K. Although the topological invariance of 4>^{h) is not known, 
from the computational point of view, this gives an excellent consistency check of correctness 
of the code. 

Equations (2a)-(2c) are manifest by the definition of In [DG13] it was shown that 

T.y is a topological invariant, defined up to a sign. We may call the 1-loop invariant. If we 
write 

( OO 

n=2 

then are the n-loop invariants of the 7 . In [DG13] it was conjectured that is well- 
dehned up to addition of an integer multiple of 1/24, and that are topological invariants 
for n > 3. 

The dehnition of is given explicitly by formal Gaussian integration. It follows that 

S-y^n is a weighted sum of a hnite set of Feynman diagrams with Feynman loop number 
at most n. The Feynman rules were explained in detail in [DG13, Sec.1.6-1. 8 ] and the 
contributing Feynman diagrams for n = 2 and n = 3 were explicitly drawn. For n > 3, the 
number of Feynman diagrams gets large and drawings-by-hand is not advisable. 

For the beneht of the reader, we recall the Feynman rules from [DG13, Sec.1.6-1. 8 ]. By 
connected Feynman diagram G we mean a connected multigraph, possibly with loops and 
multiple edges. If G is a Feynman diagram, its Feynman loop number L{G) is given by 

L(G) = |Fi(G)| + |F 2 (G)|+ 61 (G) 

where |I4(G)| is the number of fc-valent vertices of G and 61 (G) is the first betti number (also 
known as the number of holes) of G. It is easy to see that a connected Feynman diagram 
with loop number at most n has at most 2n — 2 vertices and at most n holes. Hence, there 
are hnitely many Feynman diagrams of loop number at most n. 

Fix a Neumann-Zagier datum 7 = (A,B, I/, /, /", z) which we assume is non-degenerate, 
that is the propagator (defined below) makes sense. In each Feynman diagram G, the edges 
represent an N x N propagator 

U = h (-B-^A + diag(l/(l - z)))”^ 
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n 

2 

3 

4 

5 

6 

9n 

6 

40 

331 

3700 

53758 


Table 1. The number gn of graphs that contribute to the n-loop invariant for n = 2,..., 6. 


while each fc-vertex comes with an iV-vector of factors 


-p(fc) 
i ! 


p{fc) 


ak+n-L(D) , 

-i)‘ E — 




Pj \ 

- 1-^12 — n — k 


P=0!k 


Pl 


(^i ) + 


B 


k = 1 
k>2 ’ 


where = 1 (resp., 0) if /c = 1,2 (resp., k > 3). Here is the fc-th Bernoulli number 
{Bi = — 1 / 2 , i ?2 = 1 / 6 ) and Lh( 2 ;) = ^ s-polylogarithm function 

for s a nonpositive integer. The diagram G is then evaluated by contracting the vertex factors 
r) with propagators, multiplying by a standard symmetry factor, and taking the part 
of the answer. In the end, S.y^n is the sum of evaluated diagrams, plus an additional vacuum 
contribution 


p(o) 


B^ 


n\ 


N 


^Ll 

i=l 


2-niZi + 


0 


n = 2 
n > 3 


1.3. Our code. Our goal is to give an exact computation for the n-loop invariants for n = 
1,..., 6 of a Neumann-Zagier datum of a SnapPy triangulation. Our method is implemented 
in SnapPy. We accomplished this in three steps. 

(a) We wrote a Python method generate.!eynman.diagrams .py that generates all Feyn¬ 
man diagrams that contribute to the n-loop invariant. The Feynman diagrams were 
generated by hrst generating trees, and then adding to them multiple edges or loops. 
The number of such diagrams is shown in Table 1. Observe that if G is a multi¬ 
graph with corresponding simple graph S{G) then S{G) has at most 2 n — 2 ver¬ 
tices and at most n holes, and L{G) can be obtained from S{G) by adding at most 
n — L{S{G)) + |Hi(S'(G))| -I- |V' 2 (S'(G))| edges. Thus, all Feynman diagrams with 
Feynman loop number at most n can thus be generated by hrst generating all trees 
with at most 2n — 2 vertices then iteratively adding edges between pairs of vertices. 
Every edge added also adds an additional hole. If multigraph G has more than 
n — |Vi(G)| — |V 2 (G)| holes it cannot be the subgraph of a Feynman diagram with 
Feynman loop number at most n. 

(b) We wrote a Python class NeumannZagierDatum which gives the Neumann-Zagier 
matrices and the exact value of the shape parameters that recover the geometric rep¬ 
resentation of an ideal triangulation. The exact computation of the shape parameters 
was done using the Ptolemy module [GGZ15, GDW] and the numerical computation 
is already implemented in SnapPy. 

(c) We wrote Python classes nloop.exact. py and nloop_num. py which given a Neumann- 
Zagier datum 7 and a natural number n = 1,..., 6 computes exactly (as an 
element of the trace held) or numerically to arbitrary precision. 
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To verify correctness of our code, we computed the u-loop invariants for n = 1,..., 5 for 
different triangulations of each of a hxed knot, such as 52 , (—2, 3, 7) pretzel, 6 i, and 62 . In all 
cases the results agreed (up to a sign when n = 1 and up to addition of 1/24 times an integer 
when n = 2). This illustrates both the topological invariance of the n-loop invariants, and 
the correctness of our code. 

1.4. Usage. The essence of our code lies in two Python classes NeumannZagierDatuin and 
nloop. The former takes as input a manifold and generates the Neumann-Zagier datum 
7 = [A, B,u, f, f", z), and the latter takes as input Neumann-Zagier datum, an integer n, 
and a list of Feynman diagrams and returns the n-th loop invariant S^^n- 

The NeumannZagierDatuin class has three optional arguments engine, verbose, and f ilemame 
which are set to None, False, and None, respectively, by default. The engine variable is 
passed as an option into the Ptolemy module and controls the method in which solutions 
to the Ptolemy variety are found. The preferred value for this variable for manifolds in 
CensusKnots is engine= "magma", which refers to the Sage interface to the Magma Com¬ 
putational Algebra System [BCP97]. If Magma is not available, engine="None" will at¬ 
tempt to compute solutions of the Ptolemy variety using Sage. Solutions for manifolds 
in HTLinkExteriors and LinkExteriors have been precomputed and are available with 
the Ptolemy module using engine="retrieve", [GGZ15]. This option requires an internet 
connection, but will automatically switch to recomputing locally if the download is unsuc¬ 
cessful. The output of the Ptolemy module including the retrieve option are suppressed 
with verbose=False and are displayed with verbose=True. 

To utilize the NeumannZagierDatum class use a terminal to navigate to the directory 
containing nloop.exact .py, available at [GSS], and load Sage. Once loaded, the class must 
hrst be initiated via 

sage: attach(’nloop.exact .pyO 
sage: M = Manifold('6_1 ’) 

sage: D = NeumannZagierDatum(M, engine="retrieve"). 

To generate the Neumann-Zagier datum use 
sage: D. generate_nz_data(). 

This will assign a Python list [A, B, u, /, /", z, embedding] consisting of the Neumann-Zagier 
datum plus the embedding of the something in the something to the class variable nz. If 
the optional argument file_name is used, this variable will be saved as a Sage object hie 
(*.sobj) in the current directory. To view the data simply use 
sage: D.nz. 

The shape equations z and held embedding may be computed separately via 
sage: D. exact_shapes_via_ptolemy_lifted() 
and 

sage: D. compute_ptolemy_f ield_and_embedding(), 
respectively. 

Once the Neumann-Zagier datum has been computed, one may use it to compute the 
n-loop invariants S^^n- First, load the Feynman diagrams you wish to use and choose an 
invariant you wish to calculate, 
sage: n = 2 
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sage: diagrams = load('Gdiagrams.sobj’) 
sage: E = nloop(D.nz, n, diagrams). 

Here, we have chosen to calculate using Feynman diagrams up to six loops. Note that 
the manifold M is not directly used when initiating the nloop class, as all the information 
about the manifold we need is encoded in the Neumann-Zagier datum D.nz. To compute 
the invariant use 

sage: E.one_loop() 
if n = 1 or 

sage: E.nloop_invariant() 

otherwise. To do this using a precomputed Neumann-Zagier datum Sage object hie instead 
of dehning D as above use 

sage: nz = load(^nz_exact_6_l.sobj O 
sage: E = nloop(nz, n, diagrams). 

The entire process described above has been streamlined into two automated functions for 
convenience. For example, to start with a specihed manifold M and diagrams list, compute 
the Neumann-Zagier datum, and then compute the n-loop invariant, simply use 
sage: nloop_from_manifold(M, n, diagrams, engine="retrieve"). 

The NeumannZagierDatum optional arguments described above may be entered here as seen 
in the example. On the other hand, to start with a precomputed Neumann-Zagier datum 
Sage object hie (loaded as nz) and a diagrams list then compute the n-loop invariant, simply 
use 

sage: nloop_from_nzdatum(nz, n, diagrams, engine="retrieve"). 

Also available at [GSS] is an almost identical version of our code, nloopmum.py, which 
produces numerical results to arbitrary precision instead of exact computations. The usage 
for this hie is the same. 


1.5. Sample computations. The results of our computations are available from [GSS], 
along with the code and data hies. 

To illustrate our method, consider the 6 i = K4i knot with trace held = Q(a:), where 
X = -1.50410836415074 ■ • • + il.22685163774658 ... is a root of 


-I- 2x^ -I- x^ — 3x -I- 1 = 0. 


is a number held of type [0, 2] with discriminant 257, a prime number. It follows that 
the Bloch group is an abelian group of rank 2 [SusOO, ZicOOj. The default SnapPy 

triangulation for K4i uses 4 ideal tetrahedra with shapes 


.’33720 5^30. o I332 313820 1 

'22 2’ ’22 2’ 22 2 


A Neumann-Zagier datum 7 = (A, B, v, /, /", z) is given by 



/ 1 0 -1 0 \ 


/ 1 0 0 1 \ 


( 1 ^ 


( 1 ^ 


(^\ 


-11 11 


0 10 1 


2 


2 

nil 

0 

A = 

-11 10 

; B — 

0 0 10 


1 

, / = 

0 

, / = 
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V-1 1 0 W 


\ 0 0 0 2 y 


V2y 


V1 y 


VO/ 
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The n-loop invariants for n = 1,..., 

17 


6 are given by: 


■ 3 17 2 

r = --x^ - X 

2 2 


7 
2 

46490 3 


-X + 6 


3 231209 2 

^ “ 198147^ 396294^ 792588^ ~ 264196 

570416 3 2833463 2 1122215 1386486 

^ “ 16974593^ 33949186^ 16974593^ ~ 16974593 

2255130587026 3 91695358340911 2 85651263871967 1596902056811 

^ “ ”50451970187565^ ~ 807231523001040^ ~ 807231523001040^"^ 20180788075026 
^ 37040877003091 3 330280282463219 2 53499149965837 72838757049049 

® ~ ~ 1728890845093894^ ~ 6915983380375576^ ~ 1798890845093894^”'' 115954793006959 


473191 62777 

264196 


50451970187565" 

^ 37040877003091 ^ _____ __ _ ^ __ __ _ 

® “ ”1728820845093894^ ” 6915283380375576^ ” 1728820845093894" 
1449319256564305241317 3 23592842410230239076799 

® “ 17984434859623040256945 ^ 115100383101587457644448^ 

110567432832899754708187 ^ 20008494585620168748319 

575501915507937288222240^ ” 143875478876984322055560 


1152547230062596 
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